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Abstract 

Anisotropic decompositions using representation systems based on parabolic scaling such as curve- 
lets or shearlets have recently attracted significantly increased attention due to the fact that they 
were shown to provide optimally sparse approximations of functions exhibiting singularities on lower 
dimensional embedded manifolds. The literature now contains various direct proofs of this fact and 
of related sparse approximation results. However, it seems quite cumbersome to prove such a canon 
^ ^ of results for each system separately, while many of the systems exhibit certain similarities. 

P5 In this paper, with the introduction of the notion of parabolic molecules, we aim to provide 

^ a comprehensive framework which includes customarily employed representation systems based on 

' ' parabolic scaling such as curvelets and shearlets. It is shown that pairs of parabolic molecules have 

the fundamental property to be almost orthogonal in a particular sense. This result is then applied 
to analyze parabolic molecules with respect to their ability to sparsely approximate data governed 
by anisotropic features. For this, the concept of sparsity equivalence is introduced which is shown to 
allow the identification of a large class of parabolic molecules providing the same sparse approximation 
results as curvelets and shearlets. Finally, as another application, smoothness spaces associated with 
parabolic molecules are introduced providing a general theoretical approach which even leads to novel 
results for, for instance, compactly supported shearlets. 
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0\ 1 Introduction 

\^ Recently, a paradigm shift could be observed in applied mathematics, computer science, and electrical 

engineering. The novel paradigm of sparse approximations now enables not only highly efficient encoding 
04 of functions and signals, but also provides intriguing new methodologies, for instance, for recovery of 

" _ I missing data or separation of morphologically distinct components. At about the same time, scientists 

^ began to question whether wavelets are indeed perfectly suited for image processing tasks, the main reason 

being that images are governed by edges while wavelets are isotropic objects. This mismatch becomes 
also evident when recalling that Besov spaces can be characterized by the decay of wavelet coefficient 
sequences however Besov models are clearly deficient to adequate capturing of edges. 



1.1 Geometric Multiscale Analysis 

These two fundamental observations have led to the research area of geometric multiscale analysis whose 
main goal is to develop representation systems, preferably containing different scales, which are sensitive 
to anisotropic features in functions/signals and provide sparse approximations of those. Such represen- 
tation systems are typically based on parabolic scaling, and we exemplarily mention (first and second 
generation) curvelets 8J, contourlets [T^, and shearlets [25 . Browsing through the literature, it becomes 
evident that some properties such as sparse approximation of so-called cartoon images are quite similar 
for some systems such as curvelets and shearlets, whereas other systems such as contourlets show a some- 
how different behavior. Delving more into the literature we observe that for those systems exhibiting 
similar sparsity behavior many results were proven with quite resembling proofs. One might ask: Is this 
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cumbersome close repetition of proofs really necessary? We believe that the answer is no and that a 
general framework for representation systems based on parabolic scaling does not only solve this problem 
but even more provide a fundamental understanding of such systems and allow for a categorization of 
these. 



1.2 Parabolic molecules 

The main goal of this paper is to proclaim the framework of parabolic molecules as a general concept 
encompassing in particular curvelets and both band-limited and compactly supported shearlets. The idea 
of molecules in geometric multiscale analysis dates back to the seminal work by Candes and Demanet [S] , 
in which they studied the curvelet representation of wave propagators by using what they called curvelet 
molecules. Later, Guo and Labate adopted this idea and introduced shearlet molecules in |21j . 

Both such generalization approaches however suffer from the fact that they are solely designed to 
weaken the conditions of the particular respective systems, namely curvelets and shearlets. In contrast to 
this, our philosophy is to introduce molecules, which encompass a wide class of directional representation 
systems by using parabolic scaling as a unifying concept. This is justified by the fact that all known 
systems providing optimally sparse approximation of cartoon images follow a parabolic scaling law; and 
it is strongly believed that this is necessary. In fact, our framework is general enough to, for instance, 
include all known curvelet-type as well as shearlet-type constructions to date. 

Our main result (Theorem |2.9[ ) will show that the Gramian matrix between any two systems of 
parabolic molecules satisfies a strong off-diagonal decay property and is in that sense very close to a 
diagonal matrix. This will become key to transfer the celebrated properties of curvelet systems to other 
systems based on parabolic scaling; a fact which we can summarize in the following meta-result: 

Meta-Theorem. All frame systems based on parabolic scaling (specifically curvelets and shearlets) posses 
the exact same approximation properties, whenever the generating Junctions are sufficiently smooth, as 
well as localized in space and frequency. 

This meta-theorem has been verified on a case-by-case basis for a number of different systems. In 
this paper, for the first time, a rigorous framework is provided which applies to, for instance, all known 
curvelet or shearlet constructions at once. This will be exemplarily demonstrated by the results on 



sparse approximation (Theorem 4.6) and anisotropic smoothness spaces (Theorem 4.10) which are indeed 



universally applicable to all parabolic molecules. 



1.3 Sparsity Equivalence 

Focussing on the property of sparse approximation of images governed by anisotropic features, it might 
be even more beneficial to derive a categorization of parabolic molecules according to their approximation 
behavior. We accommodate this request by introducing the notion of sparsity equivalence in Subsection 
|4.1[ which leads to equivalence classes and further to the sought classification. As a byproduct, our 
framework yields a simple derivation of the results in [201 [57] from [H]. In fact, our results provide a 
systematic way to analyze the sparse approximation of cartoon images of systems by elements of the 
class of parabolic molecules categorized by equivalence classes of sparsity equivalence. 



1.4 Contribution and Expected Impact 

Summarizing, the significance of the notion of parabolic molecules as a higher level viewpoint lies in 
the fact that it not only provides a general framework which contains various directional representation 
systems as special cases and enables a quantitative comparison of such, but it moreover allows the transfer 
of results concerning properties of such systems without repeating quite similar proofs. A few examples, 
for which this conceptually new approach is fruitful, will be presented in Section [4] including optimally 
sparse approximations of cartoon-like images. 

We therefore anticipate this novel framework to have the following impacts: 

• A thorough understanding of the ingredients of representation systems based on parabolic scaling 
which are crucial for an observed behavior such as sparse approximations of cartoon images, thereby 
also categorizing different (sparsity) behaviors. 
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• A framework within which resuhs can be directly transferred from one system to others without 
repetition of similar proofs. This will allow to establish a desired result for a system based on 
parabolic scaling by choosing, for instance, a shearlet or curvelet system best suited for the proof, 
and transfer the result subsequently to any other systems by utilizing the results in this paper. 

• An approach to design new representation systems based on parabolic scaling depending on several 
parameters whose impact on, for instance, sparse approximation behavior is then known in advance. 

1.5 Extensions 

The framework introduced in this paper and the derived results are amenable to generalizations and 
extensions. We briefly discuss a few examples. 

• Other Systems. This general framework supports the introduction of novel systems based on 
parabolic scaling fulfilling a list of desiderata designed according to a particular application. Such 
systems can now be objectively compared with other systems according to, for instance, sparse 
approximation properties. 

• Systems with Continuous Parameters. One might also ask whether a similar general framework 
for systems based on parabolic scaling with continuous parameters can be introduced. In light of 
Subsection |4.1[ this however requires a different sparsity model; one conceivable path would be to 
compare their ability to resolve wavefront sets. 

• Further Properties. In this paper, we studied the impact of our general framework on the problems 
of sparse approximation and anisotropic function spaces. This strategy can certainly be also used for 
other applications such as efficient decomposition of the Radon transform, which has been studied 
both for shearlets [IS] and curvelets [7] , as well as the analysis of geometric separation as studied 
in [T3]. 

• Weighted Norms. When aiming at transferring results such as sparse decompositions of curvilinear 
integrals |^ or sparse decompositions of the Radon transform 7J, sometimes weighted £p norms 
might need to be analyzed. 

• Higher Dimensions. We have formulated our results in the bivariate setting. However, an extension 
to arbitrary dimensions is possible using essentially the same arguments. This is especially relevant 
since by now several different curvelet and shearlet constructions exist for three-dimensional data 

[301 1211123] ■ 

1.6 Outline 

This paper is organized as follows. In Section [2] the notion of parabolic molecules is introduced. It 
is then shown in Section [3] that curvelets and both band-limited and compactly supported shearlets 
are special cases of this framework. Almost orthogonality of pairs of parabolic molecules is proven in 
Sectionjs] Finally, in Section|4j this result is utilized for two applications. First, in Subsection |4.1[ using 
the novel concept of sparsity equivalence a large class of parabolic molecules providing the same sparse 
approximation results as curvelets and shearlets is identified. Second, in Subsection |4.2[ smoothness 
spaces associated with parabolic molecules are studied. 

1.7 Notation 

We comment on the notation which we shall use in the present work. Denote by Lp{M.'^) the usual 
Lebesg ue spaces with associated norm || • ||p. For a discrete set A equipped with the counting measure we 
denote the corresponding Lebesgue space by ^p(A) or £p if A is known from the context. The associated 
norm will again be denoted || • |jp. We use the symbol (•, •) indiscriminately for the inner product on the 
Hilbert space L2(M'') as well as for the Euclidean inner product on M''. The Euclidean norm (a;,a;)^/^ 
of a vector x S M'* will be denoted by For a function / e ii(M'^) we can define the Fourier 
transform f{uj) :— J^^ /(x) exp(— 27ri(x, a;))da;. By density this definition can be extended to tempered 
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distributions /. We shall also use the notation T to denote the one-dimensional torus which can be 
identified with the half-open interval [0, 27r). Sometimes wc will use the notations (a:;)+ := max(a;, 0) 
[a;J := max{; £ Z : I < x}, and {x) := (1 + x^}^/^. Finally, we use the symbol A < B to indicate that 
A < CB with a uniform constant C. 



2 Parabolic Molecules 

All anisotropic transforms based on parabolic scaling which have appeared in the literature are indexed 

by a scale parameter describing the amount of anisotropic scaling, an angular parameter describing the 
orientation and a spacial parameter describing the location of an element. Nevertheless, these specific 
constructions are based on different principles: For curvelets the scaling is done by a dilation with respect 
to polar coordinates and the orientation is enforced by rotations. Shearlets on the other hand are based 
on affine scaling of a single generator and the directionality is generated by the action of shear matrices. 

It is the purpose of this section to introduce the concept of parabolic molecules which distills the 
essential properties out of these constructions in terms of time-frequency localization properties. As it 
will turn out, all previous constructions of curvelets and shearlets are instances of this concept, a fact 
that enables us to operate in a much more general setup than in previous work. 



2.1 Definition of Parabolic Molecules 

Let us now describe our setup. We start by defining our parameter space 

P := M+ X T X 

where a point p = (s, 9,x) € P describes a scale 2*, an orientation 9, and a location x. 

Parabolic molecules are defined as systems of functions {mx)x^A with each m\ G L2(IR^) satisfying 
some additional properties. In particular, each function m\ will be associated with a unique point in 
P, as we shall make precise below. Since we are dealing with discrete systems (frames) we would like 
to operate with discrete sampling sets contained in P. We call such sampling sets parametrizations as 
defined below. 

Definition 2.1. A parametrization consists of a pair (A, $a) where A is a discrete index set and $a is 
a mapping 

A -s- P, 
AeA h-s. {sx,Ox,xx). 

which associate with each X E A a scale sx, a direction 6x and a location xx S R^. 
There exists a canonical parametrization 

A°:={{j,l,k)GZ^ : j>0, Z = -2LiJ-i, • • • , 2LiJ-i } , 
where for A = {j,l,k) we define $°(A) := {sx,9x,xx) with sx := j, Ox ■= 11~^^'^^'k, and xx ■= 

We are now ready to define parabolic molecules. Our definition essentially says that molecules have 
frequency support in parabolic wedges associated to a certain orientation and spacial support in rectangles 
with parabolic aspect ratio. 

For this, we will use the following notion. We let Re = { ^^"^m^ ) denote the rotation 



sin(^) cos(^) 

matrix of angle 6, and := diag(a, y/a) be the anisotropic dilation matrix associated with a > 0. 

Definition 2.2. Let A he a parametrization. A family {mx)xeA is called a family of parabolic molecules 

of order {R, M, Ni, N2) if it can be written as 

mx{x) = 23^-/4a(^) {D2'xRe, {x - xx)) 

such that 

dPa^^\0\ < min (1,2"^^ + \^,\ + 2-^^/'\^2\)^ m-""^ (6)"^^ (1) 
for all \P\ < R. The implicit constants are uniform over AeA. 



4 



Remark 2.3. For convenience our definition only poses conditions on the Fourier transform ofm\. The 
number R describes the spatial localization, M the number of directional (almost) vanishing moments and 
Ni , 7V2 describe the smoothness of an element mx . We also refer to Figure [I] for an illustration of the 
approximate frequency support of a parabolic molecule. 
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Figure 1: Left: The weight function min (l, 2""^ + |Ci| + '^^'^'^\^2\)^' (1^1)^^' (6)^^' for sx = 3, M = 
3, A'^i = A*'2 — 2. Right: Approximate Frequency support of a corresponding molecule rh\ with 9\ — 7r/4. 

We pause to record the following simple estimates: In polar coordinates we have the representation 

rhx{r,ip) = 2-^"'^/'^a(-^) (2-"^r cos((y3 + 61a), 2""^/^ sin((y5 + 6*^)) exp (27ri(xA, 0) • 

Equation ([ij directly implies that in polar coordinates we have the estimate 

\mx[0\ S 2-2-^/4 min (1, 2-^^ (l + r))^^(2-^^r)"^' (2"^^/^ sin(^ + (2) 

2.2 Metric Properties of Parametrizations 

In order to proceed we need to introduce some additional (metric) properties of index sets and parametriza- 
tions. The parameter space P can be equipped with a natural notion of (pseudo) distance, see [29], which 
can be extended to a distance between indices by a pullback via a parametrization. 

Definition 2.4. Following 129}/ . we define for two indices A, /i the index distance 

a;(A,/x) :=2l^--^-l(l + 2^^od(A,M)), 

and 

d(A,/i) := \9x - 0^\^ + \xx - a;^P + |(eA,a;A - a;^)|. 
where Aq — argmin{sx, s^^) and ex — (cos(0a), sin(0A))^- 

Remark 2.5. The notation uj{\, ji) is actually a slight abuse of notation since lu is acting on P. Therefore 
it should read 

a;($A(A),$M(/^)) 

for indices X € A, 11 G M with associated parametrizations ^m- order not to overload the notation 
we stick with the shorter but slightly less accurate definition. 

Remark 2.6. We wish to mention that, in fact, real-valued curvelets or shearlets are not associated with 
an angle but with a ray, i.e., 9 and O + tt need to be identified. This is not reflected in the above definition, 
which is a slight inaccuracy. The 'correct' definition should assume that \9x\ < § G P^, the projective 
line. Therefore, it should read 

d(A,/i) := \{0x-9^}\^ + \xx-x^,f + \{9x,xx-x^)\ 



5 



with {ifi} the projection of (p onto = (— 7r/2, 7r/2]. 

However, for our results it will make no difference which definition is used. Thus we decided to employ 
Definition\2.4\ which avoids additional technicalities. 



We need to impose further conditions on an index set A in order to arrive at meaningful results. The 
following definition formalizes a crucial property, which is later on required to be satisfied by an index 
set in our results. 

Definition 2.7. An index set A with associated mapping $a is called fc-admissible if 



sup U! (A, fJ.) < oo, 



and 

sup uj (A, fj.) * < oo. 

Lemma 2.8. The canonical index set A*^ is k-admissible for all k > 2. 
Proof. We aim to prove that 

sup Lu (/i, A) < oo. (3) 
Writing = j' in the definition of uj (/i, A), we need to consider 

(^l + 2"""("')d(/i,A)) \ (4) 

According to ^ Equation (A. 2)], we have 

Y (l + 29d(/i,A))-2<22(J-9)+ (5) 

for any q. Hence, for each A; > 2, Q can be estimated by 

^2-'^l^-j'l22|j'-^'l <oo, 

which proves (|3|. □ 
2.3 Main Result 

The main result of this paper essentially states that any two systems of parabolic molecules behave in the 
same way as far as approximation properties are concerned. Specifically, we show the following theorem, 
whose proof is quite technical, wherefore we postpone it to Subsection |5.2[ 



Theorem 2.9. Let {m\)\^/^, (p^)peM be two systems of parabolic molecules of order {R, M, Ni, N2) with 
R>2N, Af > 47V - ^, Ni>2N+'^, N2 > 2N. (6) 

Then 

This result shows that the Gramian matrix between any two systems of parabolic molecules satisfies 
a strong off-diagonal decay property and is in that sense very close to a diagonal matrix. As we shall 
see in Section [4j this result has a number of immediate applications, most notably for the approximation 
properties of arbitrary frames which are systems of parabolic molecules (they turn out to be equivalent!). 

We find it particularly striking that our framework is general enough to include both curvelet-type, as 
well as shearlet-type constructions (see Section [3|. Therefore, as a consequence of Theorem |2.9[ all these 
systems satisfy the same celebrated properties of the curvelet construction given in 0. To demonstrate 
the importance of our result. Section [1] discusses selected applications of Theorem 2^ such as sparsity 
equivalence and equivalence of associated smoothness spaces. 
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3 Examples of Parabolic Molecules 



Having defined parabolic molecules in Section [2] above, it is important to examine the versatility of 
this concept. This is done in the present section. The main findings are that essentially all known 
constructions in the literature can be cast in our framework and are thus amenable to the techniques and 
results developed in this paper. 

We divide the section into two subsections. In Subsection 3.1 we study so-called curvelet-like con- 
structions. These include curvelets as defined in 4J but also other constructions, such as in [3l|29]. We 
show that all these function systems are parabolic molecules. In fact, this result should not come to 
much as a surprise: In [5 a similar concept of curvelet molecules is introduced which includes all the 
above-mentioned constructions. We also show that curvelet molecules are parabolic molecules. 

The real strength of our definition of parabolic molecules is that it includes not only curvelet-type 
constructions. In fact, we consider it one of the main findings of this paper that also shearlet-type systems 
can be thought of as instances of parabolic molecules, associated to a specific shearlet parametrization 
$'^! We show this result, as well as the admissibility of below in Subsection |3.2[ After that, to provide 
some concrete examples, we study several specific constructions. In particular, we show that compactly 
supported shearlet constructions (see e.g. [23]) are parabolic molecules. 



3.1 Curvelet-like constructions 
3.1.1 Second Generation Curvelets 

It is easily verified that curvelet molecules as defined in 5} are instances of parabolic molecules associated 
with the canonical parametrization. In particular, second generation curvelets |1] are parabolic molecules 
of arbitrary order. We start by describing the construction. Pick two window functions W{r), V{t) which 
are both real, nonnegative, C°° and supported in (5, 2) and in (—1,1) respectively. We further assume 
that these windows satisfy 

^VF(2^r)^ = l for ah reM+ and ^ V" (i - /)^ = 1 for alH G 

Now define in polar coordinates 

mo,o)(r-,w) := (2-^r) V (2^^/'iu;) and 7o-,i,fc)(-) 7(,,o,o) (• - ^u.i.k))) , 

where {j,l,k) € A". With appropriate modifications for the low-frequency case j = it is possible to 
show that the system 

r° :^{7,:AeA«} 

constitutes a Parseval frame for L2{M.'^)- In order to make the frame elements real- valued, it is possible 
to identify elements oriented in antipodal directions. This frame is customarily referred to as the tight 
frame of second generation curvelets. We now show that this frame forms a system of parabolic molecules 
of arbitrary order. 

Proposition 3.1. The second generation curvelet frame constitutes a system of parabolic molecules of 
arbitrary order associated with the canonical parametrization. 

Proof. Due to rotation invariance we may restrict ourselves to the case 9\ = 0. Therefore, denoting 
7j :— 7(j^o,o)i need to show that the function 

satisfies ([I]) for {R, AI, Ni, N2) arbitrary. First note that 

The function a^^\ together with all its derivatives has compact support in a rectangle away from the Al- 
exis. Therefore, it only remains to show that, on its support, the function a*-^-* has bounded derivatives, 
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with a bound independent of j. But this follows from elementary arguments, using r — \/ Q + Cl' 
io = arctan (^2/^1)- Then 

where 

aj(0 := 2-^Y22j^fT2^, and l3j{^) := 2^'/^ arctan 

Now it is a simple calculus exercise to show that all derivatives of aj and Pj are bounded on the support 
of a^^^ and uniformly in j. This proves the result. □ 

3.1.2 Hart Smith's Parabolic Frame 

Historically, the first instance of a decomposition into parabolic molecules can be found in Hart Smith's 
work on Fourier Integral Operators and Wave Equations [22] • This frame, as well as its dual, again forms 
a system of parabolic molecules of arbitrary order associated with the canonical parametrization. We 
refer to [291 H] for the details of the construction which is essentially identical to the curvelet construction, 
with primal and dual frame being allowed to differ. The same discussion as above for curvelets shows 
that this system consists of parabolic molecules. 

3.1.3 Borup and Nielsen's Construction 

Another very similar construction has been given in [3j. In this paper, the focus has been on the study of 
associated function spaces. Again, it is straightforward to prove that this system constitutes a system of 
parabolic molecules of arbitrary order associated with the canonical parametrization. As a corollary to 
our results, it will actually turn out that the spaces defined in ^ coincide with the approximation spaces 
corresponding to curvelets, shearlets, and Smith's transform. 

3.1.4 Curvelet Molecules 

In [S] the authors introduced the notion of curvelet molecules which are a useful concept in proving 
sparsity properties of wave propagators. For the sake of completion, we include the exact definition. 

Definition 3.2. Let A" be the canonical parametrization. A family {m\)\^t^o is caHed a family of curvelet 
molecules of regularity R if it can be written as 

mx{x) = 23^-/4a(^) (i?2=Ai?e, {x - xa)) 

such that for all < i? and each = 0, 1, 2, . . . 

\d^a^^\x)\<{xY^ (7) 

and for M = 0,1,... 

|a(^)(0|<min (1,2-^- +161+2-^^/2161) . (8) 

This definition is similar to our definition of parabolic molecules, however with two crucial differences: 
First, ([I]) allows for arbitrary rotation angles and is therefore more general. Curvelet molecules on the 
other hand are only defined for the canonical parametrization A*' (which, in contrast to our definition, 
is not sufficiently general to also cover shearlet-type systems). Second, the decay conditions analogous 
to our condition ([T]) are more restrictive in the sense that it requires infinitely many nearly vanishing 
moments. In fact, the following result holds: 

Proposition 3.3. A system of curvelet molecules of regularity R constitutes a system of parabolic 
molecules of order (cx), 00, R/2, R/2) . 

Proof. The definition of curvelet molecules as above implies that the estimate ^ also holds for all 
derivatives of a'^^\ see [Sj. Furthermore, by ([t]), all derivatives of d^'^^ can be estimated in modulus by 
which in turn can be estimated by {\^\)~^^'^ {^2)^^^^ ■ This yields the desired estimate. □ 
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3.2 Shear lets 



Shearlets were introduced in 2006 as the first directional representation system which not only satisfies 
the same celebrated properties of curvelets, but also provides a unified treatment of the continuum and 
digital setting. This key property is achieved through utilization of a shearing matrix as a means to 
parameterize orientation, which is highly adapted to the digital grid in contrast to rotation. For more 
information on shearlets, we refer to the book |25j . 

It is perhaps not surprising that curvelets and their relatives described above fall into the framework 
of parabolic molecules. Here we show the crucial fact that shearlets can be seen as a special case of 
parabolic molecules as well. Consider the discrete index set 

A'' := {{e,j,l,k) e X : e e {0,1}, j>0, I = -2^iK ■ ■ ■ ,2^^^} , (9) 

and the shearlet system 

E := {tJA : A e A^} , 

with 

1 /2-Lj72J 



where = Da, D\ := diag(v^, a), Sij •= ( q ^ ) ^^'^ ^ V'^J ' "^^^^ define shearlet 
molecules of order {R, M, Ni,N2), which is a generalization of shearlets adapted to parabolic molecules. 



in particular including the classical shearlet molecules introduced in [21], see Subsection 3.2.5 



Definition 3.4. We call S a system of shearlet molecules of order (R, M, Ni, N2) if the functions 
V, ^li,k^ i^],i,k sa-i^sfy 

\d^r,,A^i.^2)\ < niin(l, |6+e|)*' (ier^M6-e)^^^ (10) 
for every (3 € with \(3\ < R. 

Remark 3.5. In our proofs it is nowhere required that the directional parameter I runs between — 2L3J and 
—2^i^ Indeed, I running in any discrete interval —C2^i^ , ■ ■ ■ , C2L5J would yield the exact same results, as 
a careful inspection of our arguments shows. Likewise, in certain shearlet constructions, the translational 
sampling runs not through k £1? hut through t1? with t > a sampling constant. Our results are also 
valid for this case with the similar proofs. The same remark applies to all curvelet-type constructions. 

Now we can show the main result of this section, namely that shearlet systems with generators 
satisfying ( 10 ) are actually instances of parabolic molecules associated with a specific shearlet- adapted 
parametrization 

Proposition 3.6. Assume that the shearlet system E constitutes a system of shearlet molecules of order 
{R, M, Ni, N2). Then E constitutes a system of parabolic molecules of order (R, M, Ni, N2), associated 
to the parametrization (A°',<i>°'), where 

$'^(A) = {sx,9x,xx) ~ (j,£7r/2 + arctan(-;2-L:'/2j)^(^s)-i^e_^.^^ 

Proof. We confine the discussion to e = 0, the other case being the same. Further, we will suppress the 
superscript e as well as the subscript j,l,k in our notation. We need to show that 

a(^)(.) ■.^ij{D2^,Su,RlD2-.,-) 

satisfies ([l]). The Fourier transform of a*^^' is given by 

d^'\-) = i^(D2-s,SiZRlD2s,-). 



The matrix S, J Rj has the form 



TpT _ f cos{0x) sin{0x) \ _ f a b 

'■^^ «^ I -l2-^'^/^i sini0x) +cos{0x) ) '\Q c 
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We claim that the quantities a and c are uniformly bounded from above and below, independent of j, 
To see this, consider the functions 

T(a;) := cos(arctan(a;)) and := a; sin(arctan(a;)) + cos(arctan(2;)), 

which are bounded from above and below on [—1,1], as an elementary discussion shows (in fact this 
boundedness holds on any compact interval). Clearly, we have 



and 



< 1, which, by the above implies 



Since we are only considering indices with e = 0, we have —12^ 
uniform upper and lower boundedness of the quantities a, c, i.e., there exist numbers < da < < oo, 
< 5c < Ac < oo such that for all j,l we have 



The matrix D2-^!,RJ^S^J^D2^x 



5a o, < Aa and 5c < c < Ac 
has the form 



c 



Using the upper boundedness of a, &, c and the chain rule, we can estimate for any < R: 



< 



sup 

\i\<R 



c 



< 



M 



For the last estimate we utilized the moment estimate for -0, which is given by (10). This gives us the 
moment property required in ([T]). 

Now we need to show the decay of d^a^^^ for large frequencies ^ Again, due to the fact that a, 6, c 
are bounded from above and a, c from below, and utilizing the large frequency decay estimate in (10), we 
can estimate 



< 



sup 

|7|<« 



c 



< 



a 2-'>-/^b 
c 



-Ni 



(c6) 



-N2 



The statement is proven. 



□ 



The following result shows that, just like the canonical parametrization, the shearlet parametrization 
A°' is admissible. 

Lemma 3.7. The shearlet parametrization (A*^, $'^) is k-admissihle for k > 2. 

Proof. We show the analogue to Equation ^ for the shearlet parametrization, the rest of the proof is 
analogous to the proof of Lemma |2.8[ Hence, we aim to prove that 

(11) 



J2 (l + 29d(M,A))-2 <220'-9)+ 
\eA''.s^=j 



for any q and /i G A°. Without loss of generality we assume that 0^ = 0, cc^^ = (the general case follows 
identical arguments with slightly more notational effort). Further, as before we only restrict ourselves to 
the case e = 0, the other case being exactly the same. 

First we consider the case q > j. In this situation, the expression in ( 11 ) can be bounded by a uniform 
constant. 

Now we turn to the other case j > q. In this case we use the fact that, whenever \l\ < 2^-'/^, we have 



arctan (-^^^i^^^ 



> 



/2-L4J 



and \S^^D2-,k\ > \D2-,k\, 
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to estimate ( [TI| ) by 



This can be interpreted as a Ricmann sum and bounded (up to a constant) by the corresponding integral 



dx 



dy 

2-3i/2 L 2-J/2 



(l + 2«(y2+a;2 + 1x1^X22/1))" 



compare [5| Equation (A. 3)]. This integral is bounded by a constant times 2^*^^ as can be seen by 
making the substitution Xi — >■ 2'^xi, X2 — >■ 2'^^^X2, y This shows (111 and thus completes the 

proof. □ 

These results show that the parabolic molecule concept is a unification of previous systems. In the 
remainder of this section we examine the shearlet constructions which are on the market and show that 
they indeed fit into our framework. 

3.2.1 Bandlimited Shearlets 

We start with the classical shearlet construction which yields bandlimited generators. We consider two 
functions V'l, "02 satisfying 



supp Tpl C 



1 


1 ■ 


U 


■ 1 1" 


72'" 


16. 


16' 2 



El^i (2- 



1 for \uj\ > 



supp -02 C [—1,1] 
1 



and 



j>0 

2L3/2J 



/=_2Lj/2J 

Now we define our basic shearlet via 



1 for \uj\ < 1. 



01 (6)^2 

It follows from standard arguments that the system 

SO := {23J/V {D^2.Sl • -fc) : j > 0, / = -2L4J, 
constitutes a Parseval frame for the Hilbert space L2 (C)^ with 

C:^{«:|«,l>i |{<1 



In the same way we can construct a Parseval frame for L2 (C 

C':^{M6l>i.M<l 

by reversing the coordinate axes. Finally, we can consider a Parseval frame 

$ := {ip{- -k): kel?) 

I 2\ ^ 

for the Hilbert space -^2 ( [— |, |] ) • 

Proposition 3.8. The system S := S^'UE^U^ constitutes a shearlet frame which is a system of parabolic 
molecules of arbitrary order. 

Proof. To show this, by Proposition 3.6 all we need to show is that the generators ip'^, 0^ satisfy (10 1 for 
arbitrary orders {R, M, Ni, N2). This, however, follows directly from the fact that the underlying basis 
functions are bandlimited. □ 
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3.2.2 Bandlimited Shearlets with Nice Duals 

The bandlimited shearlet frame S as described above suffers from the fact that we do not know much 
about its dual frames. In particular, we do not know whether there exists a dual frame which is also a 
system of parabolic molecules. For several results such as those in Subsection |4.2| it is however necessary 
to have such a construction. In [TS] this problem was successfully resolved by carefully glueing together 
the two bandlimited frames associated with the two frequency cones. In other words, there exist shearlet 
frames S with dual frame E' such that both S and S' form systems of parabolic molecules of arbitrary 
order. 



3.2.3 Smooth Parseval Frames of Shearlets 

In 122] another modification of the bandlimited shearlet construction is given by carefully glueing together 
two boundary elements along the seamlines with angle 7r/4. It can be shown that this yields a Parseval 
frame with smooth and well-localized elements. Again, it is straightforward to check that the system 
constructed in [32] constitutes a system of parabolic molecules of arbitrary order. 



3.2.4 Compactly Supported Shearlets 

We next analyze compactly supported shearlets [211 , and prove that they also form instances of parabolic 
molecules. Currently known constructions of compactly supported shearlets involve separable generators, 
i.e., 

1p°ixi,X2) ■■= 1plixi)ij2{x2), 1p^ixi,X2) ■■= 1p°ix2,Xi). (12) 



with a wavelet ipi and a scaling function ^^2 • We would like to find conditions on tpi , 'tp2 such that ( |10[ ) is 
satisfied for given parameters {R, M, Ni, N2), i.e., that the associated shearlet frame forms a system of 
shearlet molecules. 

First we define the crucial property of vanishing moments for univariate wavelets. 
Definition 3.9. A univariate function g possesses M vanishing moments if 

g{x)x^dx = 0, for allk^Q,...,M ~ 1. 



In the frequency domain, vanishing moments are characterized by polynomial decay near zero, as is 
well known. 

Lemma 3.10. Suppose that g : M — > C is continuous, compactly supported and possesses M vanishing 
moments. Then 

l5(0l<nun(l,|e|)"^ 

Proof. First, note that, since g is continuous and compactly supported, it is in ii(M) and therefore its 
Fourier transform is bounded. This takes care of frequencies ^ with > 1. For small ^ observe that, up 
to a constant we have 

k 



g{x)x'^dx=(^^^ 5(0) 



Hence, if g possesses M vanishing moments, all derivatives of order < M of the Fourier transform g 
vanish at 0. Furthermore, since g is compactly supported, its Fourier transform is analytic. Therefore 

m)\<\i\'\ 

which proves the claim. □ 
Proposition 3.11. Assume that tpi e is a compactly supported wavelet with M + R vanishing 



moments, and ip2 G Q^igg compactly supported. Then, with "0^ defined by (12), the associated 

shearlet system E constitutes a system of parabolic molecules of order {R, M, Ni, N2). 
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Proof. In view of Proposition |3.6| we need to show that the estimate ( 10 1 holds. We only consider the 
case £ = and drop the superscript. The inverse Fourier transform of d^t/j is, up to a constant given by 
x^iIj{x). We first handle the case |^i| > 1. By smoothness and compact support of ipi,tp2, we find that 
for any \P\ < R the function 

is in Li(]R), hence it has a bounded Fourier transform which is given, up to a constant by 
It follows that the function 

is bounded. Using the simple fact that {x){y) < {\/ x^ + y"^), we get 
Now let /3 be such that |/3i| < R. Then the function 



restricted to the variable xi possesses at least M vanishing moments, due to the assumption that ip- 



possesses M + R vanishing moments. Lemma 3.10 then proves the decay of order min (l, for the 



derivatives of ip. □ 

Remark 3.12. Several assumptions on the generators ip^ , ■0^ could he weakened, for instance the sep- 
arability of the shearlet generators is not crucial for the arguments to go through. In particular, our 
arguments nowhere require neither compact support nor bandlimitedness. 

3.2.5 Shearlet Molecules of Hi] 



In [21] the results of [5] are established for shearlets instead of curvelets. A crucial tool in the proof is 
the introduction of a certain type of shearlet molecules which are similar to curvelet molecules discussed 
above, but tailored to the shearing operation rather than rotations. 

Definition 3.13. Let A°" be the shearlet index set as in A family (mA)AeA" is called a family of 
shearlet molecules of regularity R if it can be written as 

mx{x) = 23^^/'*a(-^) {D'^s^S^^x - k) 

such that for all |/3| < R and each N — Q,l,2, . . . 

\d^a'^^\x)\ < (x)^^ 

and for M = 0,1,... 

|a(^)(0| <min (l, 2"^^ + |a| + 2-^^/21^1 
By the results in |21j, the shearlet molecules defined therein satisfy the inequality ( |lO|) with the 



M 



choice of parameters {R, N, Ni, N2) ~ (00, 00, i?/2, i?/2). Therefore, in view of Proposition 3.6 shearlet 



molecules of regularity R as defined in [21 form systems of parabolic molecules of order (00, 00, R/2, R/2). 



Thus, we derive an analogous result to Lemma 3.3 for shearlet molecules: 



Proposition 3.14. A .system of shearlet molecules of regularity R constitutes a system of parabolic 
molecules of order (00, 00, i?/2, R/2) . 

Let us finish this section on examples of parabolic molecules by making the following 

Remark 3.15. By now we hope to have convinced the reader that there is a whole zoo of different 
constructions in the literature which can all be put under one roof using the concept of parabolic molecules. 
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4 Applications 



In this section we discuss selected applications of the developed theory. A particular focus will be on 
approximation properties of parabolic molecule systems {mx)x^/^, especially if they form a frame, e.g., 



AeA 

see e.g., [H]. It is well known, that in this case one can robustly represent any function / e L2{M.'^) as a 
sum 

AeA 

where {mx)x<£A is a dual frame. Approximation properties of a frame system (r7iA)AeA arc usually studied 
in terms of the sparsity of the coefficient sequence ((/, 'm'\))xeA- Below, in Subsection 4.1 we show that 



essentially any frame system which consists of parabolic molecules satisfies the same approximation prop- 
erties as the curvelet frame constructed in [3]. This, in particular, implies almost optimal approximation 
results for the class of cartoon images (see below for a definition) for all constructions mentioned in Sec- 
tion [3j for instance compactly supported or bandlimited shearlets. The above-mentioned approximation 
property may actually be regarded as the main raison d'etre of curvelet-like systems and is therefore of 
central importance. 

In Subsection |4.2| we go further and show that practically any reasonable definition of a function space 
norm based on a coefficient sequence ((/, ™a))a£a equivalent for any two frame systems consisting of 
parabolic molecules. This shows for instance that finiteness of a function space norm defined via the 
curvelet frame implies finiteness of the analogous norm defined via compactly supported shearlet frames, 
whenever the generators possess sufficient smoothness and directional vanishing moments. This result 
has not been known before. 

Remark 4.1. It is in general not the case that the dual frame {rhx)xi£A of a frame (rnx)xeA of parabolic 
molecules needs to consist of parabolic molecules, too. However, it can be shown, based on the concept of 
intrinsic localization, that the so-called canonical dual frame of {mx)x£A is of a similar form in a certain 
sense 119/ . 



4.1 Sparse Image Approximation 

Multivariate problems are typically governed by anisotropic features such as edges in images. A cus- 
tomarily employed model for such data is the class £'^{M.'^) of so-called cartoon images which is defined 

by 

where B C [0,1]^ with dB a closed C^-curve and /o,/i G Co([0, 1]^). Questions of efficient encoding 
of such a model class can be formulated in terms of optimal approximation properties. Given a frame 
system {mx)x ^ -^2(1R^), an appropriate measure for the approximation behavior is the decay rate of the 
error of best N-term approximation, i.e., of ||/ — /AfHii where /^r denotes the best iV-term approximation 
by (toa)a of some / e £^(M^), obtained as 

/at = argminll/ - ^ cxmx\\l s.t. #An < N. 

AeAjv 

A small technical problem occurs due to the fact that the representation system might not form an 
orthonormal basis in which case the computation of the best TV-term approximation is far from being 
understood. To circumvent this problem, usually the error of approximation by the N largest coefficients 
of ((/, "^A))AeA considered, which then certainly also provides a bound for the error of best A^-term 
approximation. Typically, the asymptotics of this error are studied in terms of the ^p-norms of the 
coefficient sequences of / for small values of p. Indeed, it is easily seen that membership of the coefficient 
sequence of / in an ip space for small p implies good A'^-term approximation rates, whenever the given 
representation system constitutes a frame, see, e.g., p7[ ITU] . 
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In |14| it was shown that the optimally achievable decay rate of the error of approximation of some 
/ G under the natural assumption of polynomial depth search is 

\\f - fN\\l-N~^, asiV^oo. 

Furthermore, it was proven in [5| and in [5D1 [^7] that both curvelets and shearlets attain this rate up to 
a log-factor. Apparently, these (parabolic) systems behave similarly concerning sparse approximation of 
anisotropic features. 

The next definition provides a formalization of this concept by introducing the notion of sparsity 
equivalence. It is based on the connection between best iV-term approximation rate and £p norms. 

Definition 4.2. Let {mx)xeA md {Pfj)^j,eM be systems of parabolic molecules of order (R, M, Ni, N2) 

and (R^ M , Ni, N2), respectively, and letO<p<l. Then {m\)\^A o,nd {p^)fj,eii 'f'^ sparsity equivalent 
in (.p, if 



(("^A,PM>)AeA,peAO 



< 00. 



Intuitively, systems of parabolic molecules being in the same sparsity equivalence class should have 
the same sparse approximation behavior with respect to cartoon images. The next result shows that this 
is indeed the case. 

Proposition 4.3. Let {m\)\^/^ and {Pfj)^£M be systems of parabolic molecules of order (R, M, Ni, N2) 

and {R, M , Ni, N2) , respectively, which are sparsity equivalent in £2/3- If i'm-\)\i£A possesses an almost 
best N-term approximation rate of order N^^^^ for cartoon images for any e > 0, then so does {p^)^izm- 

Proof. This is a direct consequence of the definition of sparsity equivalence and standard arguments, see 
for instance [TU]. □ 

This result enables us to provide a very general class of systems of parabolic molecules which optimally 
sparsely approximate cartoon images by using the known result for curvelets. For this, we first analyze 
when a system is sparsity equivalent to the tight frame of bandlimited curvelets. 

First, we state a simple result concerning operator norms on discrete £p spaces. 

Lemma 4.4. Let I, J be two discrete index sets, and let A : £p{I) ip{J), p > be a linear mapping 
defined by its matrix representation A = {Ai_j).^j j^j. Then we have the bound 

l/r 

\\M\eA^)^W - ( sup ^ I A, J I'', sup ^ I Ai J I 
where r :— min(l,p). 

Proof. The proof for p <1 follows easily using the fact that 

|a + 6|P < |a|P + for a, 6 e M. 

To show the case p> \ one only shows the assertion for p — 1, 00, which is trivial. The claim then follows 
by interpolation. □ 

The next theorem proves the central fact that any system of parabolic molecules of sufficiently high 
order is sparsity equivalent to the bandlimited curvelet frame from Subsection |3.1.1| 

Theorem 4.5. Assume that < p < 1, (A, $a) is a k-admissible parametrization, and F" = (7A)AeAo the 
tight frame of bandlimited curvelets. Further, assume that {m\)\^/^ is a system of molecules associated 
with A of order {R, M, Ni,N2) such that 

R>2^, M>4^-^, 7Vi>2^ + ^, iV2>2^. 
p p 4 p 4 p 



Then {m\)\^A is sparsity equivalent to F' 
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Proof. We need to show that 



i/p 



= max sup ^ |(mA,7M)r, sup ^ \{mx,j^,)\^ < oo 



mga 



AeA" 



AeAO 



AieA 



By Theorem 2.9 we have 



It follows that 



max sup ^ |(mA,7^)|^, sup ^ |(toa,7p)I 



AieA 



AeAO 



AeA" 



peA" 



< max I sup ^ a;(A,^) sup ^ a;(A, ^) j < oo 



peA 



AeA" 



ASA" 



M6A0 



due to the fc-admissibility of the parametrization of A. 



□ 



This result in combination with Proposition 4.3 now leads to the main result of this subsection 



Theorem 4.6. Assume that {m\)\^A is a system of parabolic molecules of order (R, AI, Ni, N2) such 
that 

(i) {m\)\i^\ constitutes a frame for L2{M.^) , 
(ii) A is k-admissible for all k > 2, 
(Hi) it holds that 

R>6, A/ > 12 - ^, iVi > 6 + ^, N2>6. 

Then the frame {m\)\^\ possesses an almost best N-term approximation rate of order N^^^^ , e > 
arbitrary for cartoon images. 



Proof. This follows from Proposition 4.3 Theorem 4.5 and the fact, proven in f8|, that r° provides the 
respective A'^-term approximation rate. □ 

We remark that condition (ii) holds in particular for the shearlet parametrization. Hence this result 
allows a simple derivation of the results in [501113 from 0. In fact, Theorem 4.6 provides a systematic 
way to, in particular, prove results on sparse approximation of cartoon images. 



4.2 Function Spaces 

Based on the concept of decomposition spaces introduced in [TB] , Borup and Nielsen have studied curvelet- 
like function spaces in [3]. We would like to apply our results to show that these spaces can be char- 
acterized by the transform coefficients in any frame which also forms a system of parabolic molecules. 
Consider the curvelet frame 

r" := {7,-,,fc:(j,;,fc)eA0} 
Following [3] we define for p,q,a > 1 the function spaces G^^^ given by 



3.1.1 



introduced in Subsection 
the norm 

^ i/p\ 




\{f,lj,i,kW] ■ (13) 
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This definition might seem somewhat odd, since the summation with respect to the directional parameter 
/ is done with respect to the norm. For this reason and also for some minor technical reasons, we study 
another, similar family of function spaces, namely the spaces S'^ ^ given by the norm 



k,l 




(14) 



Remark 4.7. We would like to emphazise that all the results shown in this section also hold for the spaces 
defined by {13), but with slightly more technical effort arising from the need to handle mixed Lebesgue 
spaces 121. We also remark that the function spaces defined via (14-) can be interpreted as a decomposition 
spaces of the form studied in (31 with a mixed Lebesgue space Y — £q£p (see (31 for more information). 



For technical reasons the definition in ( 14 ) forces us to work with a slightly stronger notion of admis- 



sibility than given in Definition 2.7 



Definition 4.8. An index set A with associated mapping <i>A is called strongly (fc, ^)-admissible if it is 
k-admissible and if 

^ (1 + 2«d(^, A))-*^ < 2'(^-«)+ , 

where 

:= {A e A : sx = ]} . 

Lemma 4.9. The canonical parametrization (A°,$'^) and the shearlet parametrization (A'^,$°') are both 
strongly (fc, 2)-admissible for any k > 2. 

Proof. This has already been shown earlier in ([5| for the canonical parametrization and in (111 for the 
shearlet parametrization. □ 

The aim of this section is to show the following theorem. 

Theorem 4.10. Let S = {a\ : A e A} be a frame for L2(]R^) with dual frame S = {a\ : A S A}. 

Assume further that S, S are both parabolic molecules of arbitrary order with a strongly (fc, I) admissible 
parametrization for some k, I. Then the following are equivalent norms on Sp ^: 



l/ll 



Remark 4.11. Of course it would be possible to show a quantitative version of Theorem \4-10\ in, the 
sense that E and E are only required to form a system of parabolic molecules of finite, sufficiently large 
order, depending onp,q,a. 

Before we start with the proof of Theorem |4.10| we recall the following result which is a very useful 
inequality, sometimes called the discrete Hardy inequality, see [11) . To state this result we define for a 
sequence a — (afe)^gpj the (quasi) norm 




VfcSN 



The discrete Hardy inequalities are as follows. 

Lemma 4.12. Assume that ivith A > a and r < q we have that either 

l/r 



\b,\ <2 



-Afe 



(15) 
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or 



Then we have 



l/r 



\bk\< lEK-r 

\j=k 

|b||e < ||a|| 



(16) 



Observe that defining ak ■= || ((/j 7A))AeAfc have |l/|ls°^ = |la|Uj- Armed with these useful 

facts, we may now proceed with the proof of Theorem |4.10[ 

Proof of Theorem \4-lC\ We start by fixing some notation: 

Further, we write Tj = {7^ : /i e A^} and similar for the systems S, E. Define 

A:=(r,E), A,,j (r„I]j), A:=(S,r), A,,, (S„ T,) . 

We have 

f- = (fO"A, ff=E(fJ)"A.„ f^ = (f^)"A, f/^ = 5:(ff)"A.,. 

3>0 j>Q 

Let us first assume that p > 1. Then we would like to show that 

9\ 1/9 



whenever 



/ 




E 








/ 


f / 


E 









< 00, 



AG A, 



l/p\ 



1/9 



< 00. 



For this, we obtain 



j>0 



<E||(fJ)"A, 



di H~ 62 , 



j>0 



where 



:= ^ (fj)^ Ajj and ^ (fj)^ A, 



(17) 



Next, we will prove that the inequalities (151, (16 1 are satisfied for the sequences di, e^, respectively. By 
Lemma 4.12[ this proves the desired claim. 

We start by deriving the following estimate for df. 

rf.<EII(fJ)iipiiA..IL^^., 



To further analyze ||Aij |L , , we employ Lemma 4.4 to obtain 



|A.,j||„ < max sup V | (7^, cta)!, sup V |(7m,o'a) 



(18) 



Using the fact that F and S are parabolic molecules of arbitrary order, Theorem 2.9 implies that for A'' 
arbitrary, 

\{J,,a,)\<Loi^^,X)-''. (19) 
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By ( 19 ) and the fact that the parametrization for A is strongly admissible, we can further estimate the 



first term in (181 by 



sup ^ l.(m,A)-^=2-^I^--'I sup ^ (l + 2"""(^'^)d(A.,A)) '^<2-(^-')l*-^l. 
The second term is treated similarly, and we wind up with 

for N arbitrarily large. In particular, this implies that 



(20) 



j>i \j>i / 



l/r 



with r := min(l,g), and this is (16). Similarly we can estimate 



l/r 



e. < E 2-"^^-^ II (fj) II, = 2-"^ E 2"^' II (fJ) I ^ E (2"^' II (f. 

j<i j<i \j<i 



■J ) lip 



which is (15 1. Applying Lemma 4.12 yields 





( ( 


E 




yi>o 





q\ 1/9 



q\ 1/9 



^ IE 2"M E K/'7.)r 

, J>0 \ V PGA? 



which proves one half of the desired norm equivalence. The other half (and the case of S) can be shown 
in exactly the same way. Therefore, for p > 1, the claim of the theorem is proven. 



Let us now turn to the case p < 1. For this, we need to replace the estimate (|17p with 



i>o 



< 



< 



< 



3<t 



1/p 



l/r 



EimD 



l/r 



Ell^ 

di + e 



Ell^ 



^ir iiA- ir 



where r := mm.(p,q). Now we can use (20) and proceed as above to show that the Hardy inequalities 
are satisfied for di and e^. Then, the application of Lemma [4. 12| finishes the proof. □ 

As a corollary we can consider the shearlet frame S constructed in [TS] and briefly described in 
Subsection |3.2.2| and arrive at the following theorem. 

Theorem 4.13. The curvelet frame T'^ and the shearlet frame S constructed in fT^ span the same 
approximation spaces. 
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We remark that the same conclusion holds for the frames described in Subsections 13.1.21 and 13.1.31 
Without proof we also mention that Theorem 4.10 and Theorem 4.13 also hold for the spaces G" . The 



proof is similar but slightly more technical. 

We wish to stress that in fact this result for the first time proves the meta-theorem that curvelet- 
shearlet properties are equivalent. 



and 



Remark 4.14. A similar result to Theorem 4-13 has recently been shown in \2^ . The proofs in this 
paper only apply to bandlimited constructions which present considerably less technical difficulty. 



5 Proof of Theorem 2.9 



The present sections presents the proof of our main result, namely the almost orthogonality of any two 
systems of parabolic molecules of sufficient order. Since the argument is quite involved we start by 



collecting some useful lemmata below in Subsection 5.1 before we go on to the proof of the main result, 
Theorem [2J] in Subsection [5ll 



5.1 Some Estimates 

Here we collect several estimates which will turn out useful in the proof of Theorem |2.9| The following 
lemma can be found in [17:, Appendix K.l]. 

Lemma 5.1. For N > 1 and a, a' G M_|-, we have the inequality 

(1 + a\x\y^ (1 + a'\x - y\y^ dip < max{a,a')-\l + mm{a, a')\y\)-^ . 



As a corollary we can show the next result. 
Lemma 5.2. Assume that \9\ < ^ and N > 1. Then we have for a, a' > the inequality 



T 



{l + a\sm{Lp)\y" (1 + a'|sin((^ + 6l)|) < max(a, a')"^ (1 + min(a, a^ieil)-^. (21) 

Proof. For (yS G T, we have the estimate 
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|sin((p)| > <^ lip-TT] pel2:= [f,7r] 
In order to use Lemma |5.1| we now split T into nine intervals depending on + 0, G /i, /2, /s. Then the 



left-hand side of \2l\ can be estimated by nine terms of the form 

(1 + a\^\)-'^ (1 + aV + 'J + ^l)""^ 



where d G {0, ±7r, ±27r}. By Lemma 5.1 this expression can be bounded by a constant times 

max(a, a')"^(l + min(a, a')\e + . 

Now it remains to note that for ■& G {±7r, ±27r} and \9\ < f we have |6' + i5| > \8\. This proves the 
lemma. □ 

Define the expression 

-N2 



Sx,M,N,^NAr. ^) ■■= min (l, 2"^^! + r)f' (l + 2'-/^\sm{^ + 0a) |)' 



l + 2-"V 



The following lemma will be used in order to decouple the angular and the radial variables. 
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Lemma 5.3. We have the estimate 

mill (1, + r))^' (1 + 2-^^)""^^ (l + 2-'-/\\ sin(^ + ^^A)|)"'^' < SxM-L,N,.L{r, ^) 

for every < L < N2 ■ 

Proof. After picking L we can estimate the quantity on the left hand side by 



fl 2---(l + r))^''^fl + 2--Vr'^^ r min(l,2 + \ 



L 



We need to show that 

min(l,2-*^(l + r)) 



"""^^'^ ^ < f 1 + 2-/^1 sin(^ + ^.)|r . (22) 

l + 2-''A/2r|sin(^ + 6'A)| V \ \^ ^ ^J\J \ I 



In order to prove (22), we distinguish three cases: 



r > 2^^ : In this case we derive 

min(l,2--'^(l + r)) ^ 1^1 



l + 2-''A/2r|sin(<^ + 6'A)| ~ 1 + 2--^a/2^| sin((^ + 6Ia)| ~ 1 + 2-^^/22^^ | sin(v3 + 6'a)| 

sin((/j + 9\)\ 



r < 1: For r < 1 we have 

min(l,2-'*^(l + r)) 



l + 2-^A/27,|sin((^ + 6lA)| 
• 1 < r < 2^^ : In this case we have 



< 2"^^ < 



(l + 2^^/2|sin(^ + ^^A)|) 



min(l,2-''^(l + r)) _ 1 + r 1 
l + 2--'A/2^|sin((^ + 6lA)| ^ r + 2^^/^s{n(ip + 6)\' 

Since r > 1 we have that < 2 and since r < 2^^, we have that — > 1. This proves the 
statement. 

□ 

Lemma 5.4. For A,B>0 and 

M>A-^, N2>B, Ni>A + 3/'i, 

we have the estimate 

2-i^'-+'-'> [ [ SxM,NuNAr,v)S^M,NuNAr,v)rdrd^ <2-^\'--'-\ (l + 2™"(^-^^-)/>A - ^mI) ■ 



Proof. We assume that > sx and start by showing the angular decay: By Lemma 5.2 and N2 > B, 
we have 



2-f(-A+-.) / Sx.M,N,.NAr,v)S^,M,NuNAr,v)rdrd^ < S ■ 2"^'--'-Ul + r^/^lOx - e^\) , 
where 

5:^2-2^ / min(l,2-"^(l + r))^min(l,2-''^(l + r))^(l + 2-"V)"^' (l + 2-''*'r)"^'rdr. (23) 
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The remaining estimate 



is established by sphtting up this integral into the four cases r < 1, 1 < r < 2*^, 2^^ < r < 2*'*' and 
r > 2^*^ . 

Case 1: 0<r <1 

Here we only use the moment property and estimate 



(23) < 2"^'*" 



Case 2: l<r < 2"^ 

For this case, we estimate 



-y-SxMr,-S^M2M + l, 



^ 2-s^(2+M)2-sxM 



(23) < 2- 



2 = A 



2 = A 



2-(M+2)s,,2^i''A I J. 

2"(*^+2)'*f 2^i''A2(A^+2-^i)sA 
< 2"('^+^/'*)(''^'"''^\ 



dr 



Case 3: 2"^ < r < 2"" 
For this case, we estimate 



(23) < 2 



Case 4: 2^^" < r 

For this case, we estimate 



-2s„ 



(2-^''r)''''(2-^V)"'^^rdr 



2-(2+M)sj,2^'iSA / r^-^+^~^^dr 



< 2-(2+^-f)sM2^i*A2(*-f+2-^'i)v 

_ 2^-'^i(*a'^'5a) 

<• 2-(^+3/4)(v-sa)^ 



(23) < 2 



-2s„ 



-2Wi+l 



2 — 2s,, 2-'Vi'*m2-'^i'*a . ^ 

_ 2-2s,,2^i^M2^i''^2'~^^^+^^*f 

_ 2^-'^i(^'^~*a) 

<• 2-('4+3/4)(v-''a)^ 



The proof is completed. 



□ 



5.2 Almost Orthogonality 



We now have all ingredients to prove our main result, which is Theorem 2.9 
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Proof of Theorem 2.9 To keep the notation simple, we assume that 6\ — and define sq :~ min(s>, s^). 
Further, we set 

Sx := xx- Xf^, 59 --Ox-O^. 

By definition, we can write 

where both a''^^ and b^^^' satisfy ([ij. We have the equafity 

^ 2-l(^^+^.) f (z^^^^ n^^^) Up) {D^-.^Re^i)) C'"' (exp (-27rzC • <5x)) d^,(24) 

where £ is the symmetric difi'erential operator (acting on the frequency variable) defined by 

22so Q2 



1 + 2*0 1 



12 



We have 



(exp (-27rif • Sx)) = (1 + 2^" \5x 



22so 

TT2^56'| 



{ex,5x) ) exp (— 27ri^ • (5x) . 



(25) 



for k < , we have the inequality 



where ex denotes the unit vector pointing in the direction described by the angle 9x- By Lemma 5.5 and 

R 
2 



Then, by (24) and (25) it follows that 



22so 



for aU fc < f . Now we can use Lemma 



5.4 



2'>o\Se\ 

and the fact that R > 2N to establish that 



\{mx,Pt.)\ < 2-2^1^--^-! (1 + 2^«|,50|2) ^ [ 1 + 2'°\5x\^ + 



22so 



1 + 2*0 



< 2-2^l"^-"^'l 1 + 2'<>\Se\'^ + 2'«\6x\' 



22so 



{ex,Sx) 



-N 



1 + 2*0 1 



{ex,Sx) 



< 2-^''^'^-'-^{l + 2'« {\Se\'' + \6x\^ + \{ex,6x)\)) '=uj{X,pi) 

The last inequality follows from the equation in the last line of the proof of [H Lemma 2.3]. This proves 
the desired statement. □ 

Lemma 5.5. Assume that ^ holds for two systems of parabolic molecules of order {R,M,Ni,N2). 
Utilizing the notion of the proof of Theorem \2.9[ we have 

(a(^) (Z?2-'.i?e.0^<^) {D2-^.Ro,C)) < Sxm^n,.n^MOS^,m-n,,n^^nM 

for all k < R/2. 
Proof. We show that 

-^2 



(1, 2-^-^ (1 + r))^' (1 + 2-^V) [l + 2-'-'^r\ sin(vj + Ox)]] 

•min(l,2-"^'(l + r))^'^(l + 2-"^r)"^' (l + 2-"^'/2^|sin((p + 6*^)1) (26) 
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which, usmg Lemma 5.3 with L = N2, impHes the desired statement. To show (26) we use induction 



in k, namely we show that if we have two functions a'^'''^ 6*^^^ satisfying for R, M, Ni, N2, then the 
expression 

can be written as a finite hnear combination of terms of the form 



with c, d satisfying ([T]) and R replaced by i? — 2, see Lemma [sTs] Iterating this argument we can establish 
that for k < R/2 

can be expressed as a finite linear combination of terms of the form 



(27) 



with 



< min 1,2 



+ 161 + 2-^^/^161) (iei>"^M6> 



(28) 
(29) 



and an analogous estimate for . Combining ( 28 1 and ( 29 ) , we obtain that 
11271)1 < 



min(l,2-^- + |(i52-.,i?e.0i|+2-^^/'|(i?2-=.^e.02|)''(l^2-.fifl.ei> (K^s-^.^e.OsD 
■xm^^{\,2-'- + \[D^-s,Re^i)^\+2-'-l''\[D^-s,Re^i)^^'^\\D^-.,Re^i^^^^ 

Transforming this inequality into polar coordinates as in ([2]) yields ( 26 ) . This finishes the proof. 
Lemma 5.6. Given two functions a^^^b^'^^ satisfying for R, M, Ni, Then the expression 

can be written as a finite linear combination of terms of the form 



-N2 



-N2 



□ 



with c, d satisfying ^ for R — 2, M, Ni,N2- 
Proof. Recall the definition 

To show this statement we treat the three summands of the operator £ separately. The first part is the 
identity, and therefore the statement is trivial. To handle the second part, the frequency Laplacian 2"" A, 
we use the product rule 

A(/.g) = 2 + d^^^^^ fd'^^^^^ g) + (A/).g + /(A5). 

Therefore we need to estimate the derivatives of degree 1 and the Laplacians of the two factors in the 
product 

a(^) {D2-s,Re,£,yU^^^ [D^-s.Re^O =: A{^)B{i). 
We start with the first factor, 

A{0 - a(^) (2-^^- cos(0a)6 - 2-^- sin(0A)6, 2"^^/" sin{ex)ii + 2-'-'^ cos(0a)6) • 
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Define 

:= 9(i^°)a(^) {D^-^.R9,0 and ^2(0 := d^'''^a^^'^ p2-A^e.C) • 

By definition, the functions Ai, A2 satisfy ([T]) with R replaced by R — 1. An application of the chain 
rule shows that 

= 2-^^ cos(0A)^i(e) + 2-^^/2 smi9^)A2{C)- 

Analogously, one can compute 

a(0'i)yl(^) = -2-^- sin(0A)Ai(O + 2-^-/2 cos{9^)A2{0, 
and the exact same expressions for B using the obvious definitions for Bi, i?2- We get 

5(1,0) ^^(1,0)^ = 2-'"^-"^' cos{9x)cos{9^)AiBi + 2-'^/^-"^' sm{9x) cos{9f,)A2Bi 

+ sin(0^) cos{9x)AiB2 + 2-*^/2-^m/2 sin{9x) sm{9,,)A2B2. 

It follows that 2'^° d^^'^'' Ad'^^''^'^ B can be written as a linear combination as claimed (recall that sq = 
min(sA, s^i))- The same argument applies to the product 2^'>d'--^'^'> Ad^'^'^^ B . 
It remains to consider the factor 

{AA)B + A{AB), 
where, for symmetry reasons, we only treat the summand 

{AA)B. 

In fact, it suffices to only consider 

A)B = (2-2.. cos{9xfA,i + 2-3^-/2+1 sin(0A) cos(0a)^12 - 2"^- sin(0A)'A22) B 



with Aij defined in an obvious way, satisfying (jlj) with R replaced by i? — 2. The term {d^^'^'^ A)B , and 
hence {AA)B, can be handled in the same way, as can A{AB). This takes care of the term 2'^° A in the 
definition of C. 

Finally we need to handle the last term in the definition of C, namely 

22so Q2 



l + 2«o 10^29^2 



for 9x = (otherwise the second order derivative would be in the direction of the unit vector with angle 
9\ with obvious modifications in the proof). With our notation and using the product rule we need to 
consider terms of the form 



and show that each of them, multiplied by the factor x+2^°o \9 p ' satisfies the desired representation. Let 
us start with 

which, using the fact that sin(0A) = 0, can be written as 

and which clearly satisfies the desired assertion. Now consider the expression 

which can be written as 

2-«A2-'^M cos(6'^)^iBi + 2-'>-2-'>^^^ sm{9f,)AiB2. 
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The first summand in this expression clearly causes no problems. To handle the second term we need to 
show that 

^^^^2-2-/^sin(^,)<l. (30) 
Here we have to distinguish two cases. First, assume that |(?^| < 2^-'"/^. Then we can estimate 



which readily yields the desired bound for (30). For the case > 2 we estimate 



0250 O^"" 2^"" 

-2-^^2-^/2 sin(0^) < , , , 2--°2--«/^|g^| < ^^^2-^o2-^'>/^\e^\ = 1 



which shows ( 30 ) also for this case. 

We are left with estimating the term 

which can be written as 

2-2«- cos(0^)2ABn + 2-3v/2+i sin(0^) cos(0^)^Bi2 + 2"^- sm{e^fAB2 



The first two terms are of a form already treated, and the last term can be handled using the fact that 
9,? < 01- 



sin(0^)2 < el □ 
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